Towards A Model-based Prognostics Methodology for Electrolytic 
Capacitors: A Case Study Based on Electrical Overstress 

Accelerated Aging 

Jose R. Celaya 1 , Chetan S. Kulkami 2 , Gautam Biswas 3 , and Kai Goebel 4 

1 SGT Inc. NASA Ames Research Center, Moffett Field, CA, 94035, USA 
jose. r. celaya @ nasa. gov 

2 ’ 3 Vanderbilt University, Nashville, TN, 37235, USA 
chetan. kulkarni @ Vanderbilt, edu 
biswas @ eecsmail. vuse. Vanderbilt, edu 

5 NASA Ames Research Center, Moffett Field, CA, 94035, USA 
kai. goebel @ nasa. gov 


Abstract 

A remaining useful life prediction methodology for elec- 
trolytic capacitors is presented. This methodology is based 
on the Kalman filter framework and an empirical degradation 
model. Electrolytic capacitors are used in several applications 
ranging from power supplies on critical avionics equipment 
to power drivers for electro-mechanical actuators. These de- 
vices are known for their comparatively low reliability and 
given their criticality in electronics subsystems they are a 
good candidate for component level prognostics and health 
management. Prognostics provides a way to assess remaining 
useful life of a capacitor based on its current state of health 
and its anticipated future usage and operational conditions. 
We present here also, experimental results of an accelerated 
aging test under electrical stresses. The data obtained in this 
test form the basis for a remaining life prediction algorithm 
where a model of the degradation process is suggested. This 
preliminary remaining life prediction algorithm serves as a 
demonstration of how prognostics methodologies could be 
used for electrolytic capacitors. In addition, the use degrada- 
tion progression data from accelerated aging, provides an av- 
enue for validation of applications of the Kalman filter based 
prognostics methods typically used for remaining useful life 
predictions in other applications. 
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of the Creative Commons Attribution 3.0 United States License, which per- 
mits unrestricted use, distribution, and reproduction in any medium, provided 
the original author and source are credited. 


1. Introduction 

This paper proposes the use of a model based prognostics 
approach for electrolytic capacitors. Electrolytic capacitors 
have become critical components in electronics systems in 
aeronautics and other domains. This type of capacitors are 
known for their low reliability and frequent breakdown in 
critical systems like power supplies of avionics equipment 
and electrical drivers of electro-mechanical actuators of con- 
trol surfaces. The field of prognostics for electronics com- 
ponents is concerned with the prediction of remaining useful 
life (RUL) of components and systems. In particular, it fo- 
cuses on condition-based health assessment by estimating the 
current state of health. Furthermore, it leverages the knowl- 
edge of the device physics and degradation physics to predict 
remaining useful life as a function of current state of health 
and anticipated operational and environmental conditions. 

1.1. Motivation 

The development of prognostics methodologies for the elec- 
tronics field has become more important as more electrical 
systems are being used to replace traditional systems in sev- 
eral applications in fields like aeronautics, maritime, and au- 
tomotive. The development of prognostics methods for elec- 
tronics presents several challenges due to great variety of 
components used in a system, a continuous development of 
new electronics technologies, and a general lack of under- 
standing of how electronics fail. Traditional reliability tech- 
niques in electronics tend to focus on understanding the time 
to failure for a batch of components of the same type. Just 
until recently, there has been a push to understand, in more 
depth, how a fault progresses as a function of usage, namely, 
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loading and environmental conditions. Furthermore, just un- 
til recently, it was believed that there were no precursor of 
failure indications for electronics systems. That is now un- 
derstood to be incorrect, since electronics systems, similar 
to mechanical systems, undergo a measurable wear process 
from which one can derive features that can be used to pro- 
vide early warnings to failure. These failures can be detected 
before they happen and one can potentially predict the re- 
maining useful life as a function of future usage and environ- 
mental conditions. 

Avionics systems in on-board autonomous aircraft perform 
critical functions greatly escalating the ramification of an in- 
flight malfunction (Bhatti & Ochieng, 2007; Kulkarni et al., 
2009). These systems combine physical processes, computa- 
tional hardware and software; and present unique challenges 
for fault diagnosis. A systematic analysis of these conditions 
is very important for analysis of aircraft safety and also to 
avoid catastrophic failures during flight. 

Power supplies are critical components of modern avionics 
systems. Degradations and faults of the DC-DC converter 
unit propagate to the GPS (global positioning system) and 
navigation subsystems affecting the overall operation. Ca- 
pacitors and MOSFETs (metal oxide field effect transistor) 
are the two major components, which cause degradations and 
failures in DC-DC converters (Kulkarni, Biswas, Bharadwaj, 
& Kim, 2010). Some of the more prevalent fault effects, such 
as a ripple voltage surge at the power supply output can cause 
glitches in the GPS position and velocity output, and this in 
turn, if not corrected can propagate and distort the navigation 
solution. 

Capacitors are used as filtering elements on power electronics 
systems. Electrical power drivers for motors require capaci- 
tors to filter the rail voltage for the H-bridges that provide 
bidirectional current flow to the windings of electrical mo- 
tors. These capacitors help to ensure that the heavy dynamic 
loads generated by the motors do not perturb the upstream 
power distribution system. Electrical motors are an essen- 
tial element in electro-mechanical actuators systems that are 
being used to replace hydro-mechanical actuation in control 
surfaces of future generation aircrafts. 

1.2. Previous work 

In earlier work (Kulkarni, Biswas, Koutsoukos, Goebel, & 
Celaya, 2010b), we studied the degradation of capacitors un- 
der nominal operation. There, work capacitors were used 
in a DC-DC converter and their degradation was moni- 
tored over an extended period of time. The capacitors were 
characterized every 100-120 hours of operation to capture 
degradation data for ESR and capacitance. The data col- 
lected over the period of about 4500 hours of operation were 
then mapped against an Arrhenius inspired ESR degradation 
model (Kulkarni, Biswas, Koutsoukos, Goebel, & Celaya, 


2010a). 

In following experimental work, we studied accelerated 
degradation in capacitors (Kulkarni, Biswas, Koutsoukos, 
Celaya, & Goebel, 2010). In that experiment the capaci- 
tors were subjected to high charging/discharging cycles at a 
constant frequency and their degradation progress was mon- 
itored. A preliminary approach to remaining useful life pre- 
diction of electrolytic capacitors was presented in (Celaya et 
al., 2011b). This paper here builds upon the work presented 
in the preliminary remaining useful life prediction in (Celaya 
et al., 2011a) and experimental studies done in (Celaya et al., 
2012 ). 

1.3. Other related work and current art in capacitor prog- 
nostics 

The output filter capacitor has been identified as one of the el- 
ements of a switched mode power supply that fails more fre- 
quently and has a critical impact on performance (Goodman 
et al., 2007; Judkins et al., 2007; Orsagh et al., 2005). A prog- 
nostics and health management approach for power supplies 
of avionics systems is presented in (Orsagh et al., 2005). Re- 
sults from accelerated aging of the complete supply were pre- 
sented and discussed in terms of output capacitor and power 
MOSFET failures; but there is no modeling of the degrada- 
tion process or RUL prediction for the power supply. Other 
approaches for prognostics for switched mode power supplies 
are presented in Goodman et al. (2007) and Judkins et al. 
(2007). The output ripple voltage and leakage current are 
presented as a function of time and degradation of the capac- 
itor, but no details were presented regarding the modeling of 
the degradation process and there were no technical details 
on fault detection and RUL prediction algorithms. 

A health management approach for multilayer ceramic capac- 
itors is presented in Nie et al. (2007). This approach focuses 
on the temperature-humidity bias accelerated test to replicate 
failures. A method based on Mahalanobis distance is used 
to detect abnormalities in the test data; there is no prediction 
of RUL. A data driven prognostics algorithm for multilayer 
ceramic capacitors is presented in Gu et al. (2008). This 
method uses data from accelerated aging test to detect poten- 
tial failures and to make an estimation of time of failure. 

2. Prognostics Methodology 

The process followed in the proposed prognostics method- 
ology is presented in the block diagram in Figure 1. It 
is based on a model-based prognostics framework using an 
time-dependent empirical degradation model build from ac- 
celerated aging tests. 

Accelerated Aging: The methodology is based on results 
from an accelerated life test on real electrolytic capacitors. 
This test applies electrical overstress to commercial, off the 
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shelf capacitors, in order to observe and record the degra- 
dation process and identify performance conditions in the 
neighborhood of the failure criteria in a considerably reduced 
time frame. A total of 6 accelerated aging test devices are 
available for the development of the proposed methodology. 
Electrochemical-impedance spectroscopy (EIS) is used peri- 
odically during the accelerated aging test to characterize the 
frequency response of the capacitor’s impedance. Several 
measurements are available through the aging time, includ- 
ing measurements at pristine condition and measurements af- 
ter failure condition. 

System Identification: A lumped-parameter model (M\) 
of the non-ideal capacitor impedance is assumed. This 
impedance model includes a capacitance element and an 
equivalent series resistance (ESR) parasitic element. The EIS 
measurements along with the impedance model structure are 
used in a systems identification setting to estimate the model 
parameters available throughout the aging test. This results in 
time-dependent capacitance and ESR measurements trajecto- 
ries reflecting capacitor degradation. 



Test 

Trajectory 


Health State RUL 

Estimation Estimation 


Degradation Modeling: We present here an empirical degra- 
dation model that is based on the observed degradation pro- 
cess during the accelerated life test. The objective of the 
model is to generate a parametrized model of the time- 
dependent capacitance degradation as generated by the sys- 
tem identification step. A similar degradation model can be 
generated for ESR but not considered in this work. 

Parameter Estimation: The parameters of the degradation 
model are estimated using nonlinear least-squares regression. 
The quality of the fit is good enough as to assume these pa- 
rameters as static during the prognostics process. 

Prognostics: A Bayesian framework is employed to estimate 
(track) the state of health of the capacitor based on measure- 
ment updates of key capacitor parameters. The Kalman filter 
algorithm is used to track the state of health and the degra- 
dation model is used to make predictions of remaining useful 
life once no further measurements are available. 

3. Accelerated aging experiments 

Accelerated life test methods are often used in prognostics 
research as a way to assess the effects of the degradation pro- 
cess through time. It also allows for the identification and 
study of different failure mechanisms and their relationships 
with different observable signals and parameters. In the fol- 
lowing section we present the accelerated aging methodology 
and an analysis of the degradation pattern induced by the ag- 
ing. The work presented here is based on an accelerated elec- 
trical overstress. In the following subsections, we first present 
a brief description of the aging setup followed by an analysis 
of the observed degradation. The precursor to failure is also 
identified along with the physical processes that contribute to 


Figure 1. Methodology for capacitor prognostics. 

the degradation. 

3.1. Experimental Setup 

Since the objective of this experiment is studying the effects 
of high voltage on degradation of the capacitors, the capaci- 
tors were subjected to high voltage stress through an external 
supply source using a specially developed hardware. The ca- 
pacitors are stressed under high voltage conditions and spe- 
cially developed hardware. The voltage overstress is applied 
to the capacitors as a square waveform in order to subject the 
capacitor to continuous charge and discharge cycles. 

At the beginning of the accelerated aging, the capacitors 
charge and discharge simultaneously; as time progresses and 
the capacitors degrade, the charge and discharge times vary 
for each capacitor. Even though all the capacitors under test 
are subjected to similar operating conditions, their ESR and 
capacitance values change differently. We therefore moni- 
tor charging and discharging of each capacitor under test and 
measure the input and output voltages of the capacitor. Fig- 
ure 2 shows the block diagram for the electrical overstress 
experiment. Additional details on the accelerated aging sys- 
tem are presented in (Kulkami, Biswas, Koutsoukos, Celaya, 
& Goebel, 2010). 

For this experiment six capacitors in a set were considered 
for the EOS experimental setup. Electrolytic capacitors of 
2200/iF capacitance, with a maximum rated voltage of 10V^, 
maximum current rating of 1 A and maximum operating tem- 
perature of 105 ° C were used for the study. These were 
the recommended capacitors by the manufacturer for DC-DC 
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Figure 2. Block diagram of the experimental setup. 


converters. The electrolytic capacitors under test were char- 
acterized in detail before the start of the experiment at room 
temperature. 

The measurements were recorded every 8-10 hours of the to- 
tal 180 plus hours of accelerated aging time to capture the 
rapid degradation phenomenon in the ESR and capacitance 
values. The ambient temperature for the experiment was 
controlled and kept at 25° C. During each measurement the 
voltage source was shut down, capacitors were discharged 
completely and then the characterization procedure was car- 
ried out. This was done for all the six capacitors under test. 
For further details regarding the aging experiment results and 
analysis of the measured data refer to (Kulkami, Biswas, 
Koutsoukos, Celaya, & Goebel, 2010; Celaya et al., 2011b). 

3.2. Physical interpretation of the degradation process 

There are several factors that cause electrolytic capacitors to 
fail. Continued degradation, i.e., gradual loss of functionality 
over a period of time results in the failure of the component. 
Complete loss of function is termed a catastrophic failure. 
Typically, this results in a short or open circuit in the capac- 
itor. For capacitors, degradation results in a gradual increase 
in the equivalent series resistance and decrease in capacitance 
over time. 

In this work, we study the degradation of electrolytic capac- 
itors operating under high electrical stress, i.e., V app u e d > 
Vrated • During the charging/discharging process the capaci- 
tors degrade over the period of time. A study of the literature 
indicated that the degradation could be primarily attributed 
to electrolyte evaporation, leakage current and increase in in- 
ternal pressure due to gas released due to chemical reactions 
(IEC, 2007-03; MIL-C-62F, 2008; Kulkarni, Biswas, Kout- 
soukos, Goebel, & Celaya, 2010a). An ideal capacitor would 
offer no resistance to the flow of current at its leads. However, 
the electrolyte that fills the space between the plates and the 
electrodes produces a small equivalent internal series resis- 
tance. Fig. 3 shows the structure of an electrolytic capacitor 
in detail. The ESR dissipates some of the stored energy in the 
capacitor leading to increase in the internal temperature and 
thus causing electrolyte evaporation. 

ESR and capacitance are the two main failure precursors that 
tipify the current health state of the device. ESR and capac- 
itance values were calculated after characterizing the capac- 



Figure 3. Electrolytic capacitor structure. 


itors at regular intervals. As the devices degrade due to dif- 
ferent failure mechanisms we can observe a decrease in the 
capacitance and an increase in the ESR. 

The literature on capacitor degradation shows a direct rela- 
tionship between electrolyte decrease with increase in ESR 
and decrease in capacitance value of the capacitor (Kulkarni, 
Biswas, Koutsoukos, Goebel, & Celaya, 2010b). ESR in- 
crease implies greater dissipation, and, therefore, a slow de- 
crease in the average output voltage at the capacitor leads. 

ESR and capacitance values are estimated by using a sys- 
tem identification using a lump parameter model consistent 
of the capacitance and the ESR in series as shown in Fig- 
ure 4. The frequency response of the capacitor impedance 
(measured with electro-impedance spectroscopy) is used for 
the parameter estimation. It should be noted that the lumped- 
parameter model used to estimate ESR and capacitance, is 
not the model to be used in the prognostics algorithm; it only 
allows us to estimate parameters which provide indications 
of the degradation process through time. Parameters such as 
ESR and capacitance are challenging to estimate from the in- 
situ measurements of voltage and current through the accel- 
erated aging test. 

Ci Cr Re 

0 — 1("~° — * 0 — |(“a/WW — 0 

Ideal Non ideal Capacitor 

Capacitor with parasitic series resistance 

Figure 4. Lumped parameter model (M\) for a real capacitor. 


3.3. System identification for real capacitor model 

The ESR and capacitance values were estimated from the 
capacitor impedance frequency response measured using an 
SP-150 Biologic SAS electro-impedance spectroscopy instru- 
ment. A lumped parameter model consisting of a capacitor 
with a resistor in series was assumed to estimate the ESR and 
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capacitance. 

The ideal capacitor has complex impedance Zi = 1/sCi 
where Cj is the ideal capacitance value. The complex 
impedance of model M»\ is given by 

% = R e + — — , (1) 

sC R 

where Re is the equivalent series resistance and Cr is the 
real capacitance. 

Electrochemical impedance spectroscopy measurements are 
available to characterize the electrical performance of the ca- 
pacitor. Figure 5 shows Nyquist plots of the impedance mea- 
surements for capacitor #1 at pristine condition and after ac- 
celerated aging at intervals of 71, 161 and 194 hours. The 
degradation can be observed as the Nyquist plot shifts to the 
right as a function of aging time due to increase in Re. These 
measurement are then used to estimate the parameters of the 
impedance model A4 \ from eq. (1). The parameter estimation 
performed using the EIS instrument software (EC lab). This 
is basically and optimization problem using an aggregate of 
mean squared error as an objective function. The error is ag- 
gregated at different frequencies for which measurements are 
available. The optimization is set up to minimize the objec- 
tive function by finding optimal values for Cr * and Re* • 



Figure 5. Electroimpedance measurements at different aging 
times. 

This parameter estimation is performed every time an EIS 
measurement is taken resulting on values of Cr and Re at 
different points in time through the aging of the components. 
The average pristine condition ESR value was measured to 
be 0.056 mRt and average capacitance of 2123 /jF individu- 
ally for the set of capacitors under test. 

Figure 6 shows percentage increase in the ESR value for all 
the six capacitors under test over the period of time. This 
value of ESR is calculated from the impedance measurements 


after characterizing the capacitors. Similarly, figure 7 shows 
the percentage decrease in the value of the capacitance as the 
capacitor degrades over the period under EOS test condition 
discussed. As per standards MIL-C-62F (2008), a capacitor is 
considered unhealthy if under electrical operation its ESR in- 
creases by 280 — 300% of its initial value or the capacitance 
decreases by 20% below its pristine condition value. From 
the plots in Figure 6 we observe that for the time for which the 
experiments were conducted the average ESR value increased 
by 54% — 55% while over the same period of time, the aver- 
age capacitance decreased by more than 20% (the threshold 
mark for a healthy capacitor) (see Figure 7). As a result, the 
percentage capacitance loss is selected as a precursor of fail- 
ure variable to be used in the degradation model development 
presented next. 



Figure 6. Degradation of capacitor performance, percentage 
ESR increase as a function of aging time. 



Figure 7. Degradation of capacitor performance, percentage 
capacitance loss as a function of aging time. 
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4. Degradation Modeling for Prognostics 

This sections presents the details of the degradation model 
development. A degradation model is an essential part of a 
model-based prognostics algorithm and it is typically appli- 
cation dependent. A model is formulated based on the empir- 
ical evidence of the degradation process time evolution from 
experiments presented in the previous section, particularly, 
capacitance loss as described by figure 7. 

4.1. Nominal model 


loss in capacitance. Let Ci be the percentage loss of capaci- 
tance due to degradation as shown by figure 7. The following 
equation is a degradation model V of the capacitance param- 
eter in the real capacitor model Mi . 

V x :C,(t) = e at + p, (3) 

where a and f3 are degradation model parameters that will 
be estimated from the experimental data of accelerated aging 
experiments. 


The non-ideal capacitor model Mi can be used as part of 
electronics circuits that make use of capacitors. An example 
is the low-pass filter implementation in figure 8. In this cir- 
cuit, input voltage Vi is considered as the voltage to be filtered 
and the voltage across the capacitor (this includes Re as well) 
is the output voltage V Q which is filtered. Let v(t) = V Q (t) 
and u(t) = Vi(t) in the low-pass system circuit with non- 
ideal capacitor shown in figure 8. A state- space realization 
(M 2 ) of the dynamic system is given by 


z(t) = 


C r {R + R e ) Z+ C r {R + R e ) 


u(t), (2) 


v{t) 



Re 

R + Re 


Re 


R + Re 


u(t ), 


where z(t) = Vc(t) is the state variable representing the 
capacitor voltage, Cr, Re and R are system parameters. 
Furthermore, Cr and Re are parameters that will change 
through time as the capacitor degrades. 


O — 

Vi(t) 

o — 



Figure 8. Low pass filter model (M2)- 


Model M\ describes the nominal dynamics of a low-pass fil- 
ter with a non-ideal capacitor. This model by itself is not 
sufficient to implement a model-based prognostics algorithm 
since the degradation process as reflected on model parame- 
ters is not modeled. Degradation models describing the time 
evolution of Re or Cr are needed in order to enhance Mi 
for model-based prognostics. Nevertheless, Adi is useful in 
this form for model-based fault detection and isolation which 
is not covered in this work. 


4.2. Degradation model 

The percentage loss in capacitance is used as a precursor of 
failure variable and it is used to build a model of the degrada- 
tion process. This model relates aging time to the percentage 


In order to estimate the model parameters, five capacitors are 
used for estimation, and the remaining capacitor is used to 
test the prognostics algorithm. This results in six leave-one- 
out test cases for validation of the prognostics algorithm re- 
sults. A nonlinear least-squares regression algorithm is used 
to estimate the model parameters. Table 1 presents definition 
of the test cases and the parameter estimation results. The 
estimate and 95% confidence interval is presented for param- 
eters a and (3. In addition, the error variance is included as a 
way to assess the quality of the fit. 

Figure 9 shows the estimation results for test case Tq. The ex- 
perimental data are presented together with results from the 
exponential fit function. It can be observed from the residuals 
that the estimation error increases with time. This is to be ex- 
pected since the last data point measured for all the capacitors 
fall slightly off the concave exponential model. 



Aging time (hr) 


Figure 9. Estimation results for the empirical degradation 
model. 


It should be noted that this degradation model with static pa- 
rameters will be used in a Bayesian tracking framework. This 
will help to overcome the degradation model limitation to rep- 
resent the behavior close to the failure threshold given the 
tracking framework ability to compensate the estimation as 
measurements become available. 
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Validation 

test 

Test 

capacitor 

Training 

capacitor 

a 

(95% CI) 

P 

(95% CI) 


t 2 

#2 

#l,#3-#6 

0.0162 

(0.0160, 0.0164) 

-0.8398 

(-1.1373,-0.5423) 

1.8778 

t 3 

#3 

#1, #2, #4-#6 

0.0162 

(0.0160, 0.0164) 

-0.8287 

(-1.1211,-0.5363) 

1.9654 

t 4 

#4 

#l-#3, #5, #6 

0.0161 

(0.0159, 0.0162) 

-0.8217 

(-1.1125,-0.5308) 

1.8860 

T-o 

#5 

#l-#4, #6 

0.0162 

(0.0161, 0.0164) 

-0.7847 

(-1.1134, -0.4560) 

2.1041 

T 6 

#6 

#l-#5 

0.0169 

(0.0167, 0.0170) 

-1.0049 

(-1.2646, -0.7453) 

2.9812 


Table 1. Degradation model parameter estimation results. 


4.3. State-space realization for tracking 

The estimated degradation model is used as part of a Bayesian 
tracking framework to be implemented using the Kalman fil- 
ter technique. This method requires a state-space dynamic 
model relating the degradation level at time tk to the degrada- 
tion level at time tj. c _i. The procedure to obtain a state- space 
model for T>\ is as follows. The non-linear exponential be- 
havior described in the model is represented as a first order 
differential equation which can represent the time evolution 
of Ci(t). Then, the model is discretized in time in order to 
obtain a discrete-time state-space model V 2 . 

From equation (3) we have that Ci(t ) = e at + /3, taking the 
first derivative with respect to time and substituting e at = 
Ci(t ) — /? from eq. (3) we have 

Ci = d ^f i = ac t (t) - a/3. ( 4 ) 


life of the capacitor based on the accelerated degradation data 
from the experiments described in previous sections. 

The methodology consists of the following three main steps 
and it is depicted in fig. 10. 

1. State tracking (Kalman Filter): The capacitance loss Ci 
is defined as the state variable to be estimated and the 
degradation model is expressed as a discrete time dy- 
namic model in order to estimate capacitance loss as new 
measurements become available. Direct measurements 
of the capacitance are assumed for the filter. 

2. Health state forecasting: It is necessary to forecast the 
state variable once there are no more measurements 
available at time or RUL prediction t p . This is done by 
evaluating the degradation model through time using the 
state estimate at time t p as initial value. 


Taking the finite difference approximation for Ci with time 
interval At we have 


Q(t) - Q(t - At) 
At 


aCi(t — At) — a/3, and 


Ci(t) = (1 + aA t )Ci(t - At) - af3A t . 


3. Remaining life computation: RUL is computed as the 
time between time of prediction t p and the time at which 
the forecasted state crosses the failure threshold value. 

This process is repeated for different values of t p through the 
life of the component under consideration. 


Letting tk = t and tk - i = t — At we get the state-space 
model 

Ci{tk) = (1 + aAk)Ci(tk~i) - a/3 A k . (5) 

This model can be used in a Bayesian tracking framework in 
order to continuously estimate the value of the loss in capac- 
itance through time as measurement become available. 

5. Model-based Prognostics Framework 

A model-based prognostics algorithm based on Kalman filter 
and a physics inspired empirical degradation model is pre- 
sented. This algorithm is able to predict remaining useful 





RUL(t p ) 


Figure 10. Model-based prognostics methodology 
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5.1. Kalman filter for state estimation 

A state-space dynamic model is needed for the filtering. The 
state variable Xk at time tk is defined as the percentage ca- 
pacitance loss Ci(k). Since the system measurements are 
percentage loss in capacitance as well, the output equation 
is given by yk = hxk, where the value of h is equal to one. 
The following system structure is used in the implementation 
of the filtering and the prediction using the Kalman filter. 

x k = A k x k _ 1 + B k u + V , 

, , ( 6 ) 

Vk = hx k + w, 

where, 

Ak = (1 + A k), 

B k = -a/3A k , 

h= 1, U 


The time increment between measurements is not con- 
stant since measurements were taken at non-uniform sam- 
pling rate. This implies that some of the parameters of the 
model in equation (6) will change through time. Furthermore, 
v and w are normal random variables with zero mean and Q 
and R variance respectively. The description of the Kalman 
filtering algorithm is omitted from this article. A thorough 
description of the algorithm can be found in Stengel (1994), 
a description of how the algorithm is used for forecasting can 
be found in Chatfield (2003) and an example of its usage for 
prognostics can be found in (Saha et al., 2009). 

5.2. Future state forecasting 

The use of the Kalman filter as a RUL forecasting algorithm 
requires the evolution of the state without updating the error 
covariance matrix and the posterior of the state vector. The n 
step ahead forecasting equation for the Kalman filter is given 
below. The last update is done at the time of the last measure- 
ment ti . 

n— 1 

^z+n = A n xi + A 1 B (8) 

i = 0 

The subscripts from parameters A and B are omitted since a 
constant A t is used in the forecasting mode (one prediction 
every hour). 

5.3. Noise models 

The model noise variance Q was estimated from the model 
regression residuals for each test case as presented in table 1 . 
This variance was used for the model noise in the Kalman 
filter implementation. The measurement noise variance R is 
also required in the filter implementation. This variance was 
computed from the direct measurements of the capacitance 
with the EIS equipment, the observed variance is 4.99 x 10 -7 . 


6. Prediction of Remaining Useful Life Results 

Estate estimation and RUL prediction results are discussed 
for test case T6. Figure 1 1 shows the result of the filter track- 
ing the complete degradation signal. The residuals show an 
increased error with aging time. This is to be expected given 
the results observed from the model estimation process. 



Aging time (hr) 


Figure 1 1 . Tracking results for the Kalman filter implemen- 
tation applied to test capacitor (capacitor #6). 

Figure 12 presents results from the remaining useful life pre- 
diction algorithm at time t p = 161 (hr), which is the time 
at which ESR and C measurements are taken. The failure 
threshold is considered to be a crisp value of 20% decrease in 
capacitance. End of life (EOL) is defined as the time at which 
the forecasted percentage capacity loss trajectory crosses the 
EOL threshold. Therefore, RUL is EOL minus 161 hours. 



Aging time (hr) 


Figure 12. Remaining useful life prediction at time 149 (hr). 

Figure 13 presents the capacitance loss estimation and EOL 
prediction at different points during the aging time. Predic- 
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tions are made after each point in which measurements are 
available. It can be observed that the predictions become bet- 
ter as the prediction is made closer to the actual EOL. This is 
possible because the estimation process has more information 
to update the estimates as it nears EOL. Figure 14 presents a 
zoomed-in version of figure 13 focusing in the area close to 
the failure threshold. 



Figure 13. Tq\ Health state estimation and forecasting of ca- 
pacitance loss (%) at different times t p during the aging time; 
t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 


9 












International Journal ol Prognostics and Health Management 



140 150 160 170 180 190 



140 150 160 170 180 190 



140 150 160 170 180 190 



140 150 160 170 180 190 

Aging time (hr) 


Figure 14. Tq\ Detail of the health state estimation and fore- 
casting of capacitance loss (%) at different times t p during the 
aging time; t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 


An a-X prognostics performance metric is presented in Fig- 
ure 15 for validation test T 6 . The blue line represents ground 
truth and the shaded region is corresponding to a 30% (a = 
0.3) error bound in the RUL prediction. This metric specifies 
that the prediction is within the error bound halfway between 
first prediction and EOL (A = 0.5). In addition, this metric 
allows us to visualize how the RUL prediction performance 
changes as data closer to EOL becomes available. Appendix 
B presents the a-X metric plots for the remaining validation 
cases. 



Time 


Figure 15. Performance based on a-X performance metric. 

6.1. Validation tests 

Table 2 summarizes results for the remaining life prediction 
at all points in time where measurements are available. The 
last column indicates the RUL prediction error. The magni- 
tude of the error decreases as the prediction time gets closer 
to EOL. The decrease is not monotonic which is to be ex- 
pected when using a tracking framework to estimate health 
state because the last point of estimation is used to start the 
forecasting process. 

Table 3 shows performance based on the relative accuracy 
(RA) metric in equation (9). These metrics allows for an as- 
sessment of the percentage accuracy relative to the ground- 
truth value. RA values of 100 represent perfect accuracy. The 
RA is presented for all the test cases for different prediction 
times. The last column of the Table 3 represents the me- 
dian RA of all the test cases for a particular prediction time. 
It is observed that the RA values decrease considerably for 
t p = 171. This is consistent with previous observations indi- 
cating that the algorithm with a fixed-parameter model is not 
able to cope with the sudden jump in exponential behavior 
present around the 171 hour. This is a limitation that could 
be overcome by either an enhanced degradation model or a 
an online estimation of degradation model parameters using 
a more sophisticated Bayesian tracking method like extended 
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tp 

RUL * 

RU L T2 

RU L T3 

RU L T4 

RU L T5 

RU L T q 

24 

151.04 

158.84 

164.88 

158.76 

167.76 

159.89 

47 

128.04 

131.32 

134.08 

128.35 

135.32 

125.91 

71 

104.04 

117.01 

119.88 

115.37 

122.63 

116.41 

94 

81.04 

92.69 

96.64 

93.09 

97.6 

95.42 

116 

59.04 

67.28 

65.39 

67.77 

69.5 

65.71 

139 

36.04 

44.01 

44.72 

46.88 

49.4 

53.75 

149 

26.04 

30.67 

32.41 

33.55 

35.92 

39.95 

161 

14.04 

17.23 

18.28 

18.2 

22.64 

25.6 

171 

4.04 

1.07 

2.89 

N/A 

5.52 

8.45 


Table 2. Summary of RUL forecasting results. 


Kalman filter or particle filter. 


RA = 100 



RUL * - RUL 
RUL* 


(9) 


tp 

So 

to 

ra T3 

RAt 4 

RAt5 

RAt6 

RA 

24 

94.8 

95.5 

91.9 

96.9 

99.7 

95.5 

47 

97.4 

99.3 

96.4 

96.7 

91.7 

96.7 

71 

87.5 

91.9 

84.5 

94.1 

97.1 

91.9 

94 

85.6 

90 

78.9 

94.8 

94.2 

90 

116 

86 

99.1 

76.5 

98 

96.2 

96.2 

139 

77.8 

95.8 

53.1 

96.7 

81.1 

81.1 

149 

82.1 

98.4 

46.9 

94.8 

86.6 

86.6 

161 

77.2 

87.3 

16.6 

87.5 

89.8 

87.3 

171 

26.6 

26.4 

N/A 

34.8 

63.7 

30.7 


Table 3. Validation based on relative accuracy metric. 


7. Conclusion 

This paper presents a RUL prediction algorithm based on ac- 
celerated life test data and an empirical degradation model. 
The main contributions of this work are: a) the identification 
of the lumped-parameter model (Figure 4) for a real capaci- 
tor as a viable reduced-order model for prognostics-algorithm 
development; b) the identification of the ESR and C model 
parameters as precursor of failure features; c) the develop- 
ment of an empirical degradation model based on accelerated 
life test data which accounts for shifts in capacitance as a 
function of time; d) the implementation of a Bayesian based 
health state tracking and remaining useful life prediction al- 
gorithm based on the Kalman filtering framework. One major 
contribution of this work is the prediction of remaining useful 
life for capacitors as new measurements become available. 

This capability increases the technology readiness level of 
prognostics applied to electrolytic capacitors. The results pre- 
sented here are based on accelerated life test data and on the 
accelerated life timescale. Further research will focus on de- 


velopment of functional mappings that will translate the ac- 
celerated life timescale into real usage conditions time-scale, 
where the degradation process dynamics will be slower, and 
subject to several types of stresses. The performance of the 
proposed exponential-based degradation model is satisfactory 
for this study based on the quality of the model fit to the ex- 
perimental data and the RUL prediction performance as com- 
pared to ground truth. As part of future work we will also 
focus on the exploration of additional models based on the 
physics of the degradation process and larger sample size for 
aged devices. Additional experiments are currently underway 
to increase the number of test samples. This will greatly en- 
hance the quality of the model, and guide the exploration of 
additional degradation-models, where the loading conditions 
and the environmental conditions are also accounted for to- 
wards degradation dynamics. 

Acknowledgment 

This work was funded by the NASA Aviation Safety Pro- 
gram, SSAT project. 

Nomenclature 

Cj Ideal capacitance value for an ideal capacitor 
Cr Real capacitor value for a non-ideal capacitor 
model 

Re Equivalent series resistance of the capacitor 
Ci(k) Capacitance percentage loss at time tk 
Ti Validation test on capacitor i 

Mi Nominal model for a component or system 
Vi Degradation model for a capacitor 

Rl Load resistance on electrical overstress system 
Vl Load voltage on electrical overstress system 
V 0 Electrical overstress voltage in aging system 

Zi Ideal capacitor impedance 

Z Capacitor impedance for non-ideal capacitor 

model M\ 
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A Prognostics validation results 



0 50 100 150 200 

20 





0 50 100 150 200 



0 50 100 150 200 


Aging time (hr) 

Figure 16. T 2 : Health state estimation and forecasting of ca- 
pacitance loss (%) at different times t p during the aging time; 
t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 
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Figure 17. T 2 : Detail of the health state estimation and fore- 
casting of capacitance loss (%) at different times t p during the 
aging time; t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 
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* Measured o Filtered Predicted 


^ Measured o Filtered Predicted 



Figure 18. X 3 : Health state estimation and forecasting of ca- 
pacitance loss (%) at different times t p during the aging time; 
t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 


Aging time (hr) 

Figure 19. X 3 : Detail of the health state estimation and fore- 
casting of capacitance loss (%) at different times t p during the 
aging time; t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 
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* Measured o Filtered Predicted 


^ Measured o Filtered Predicted 



Figure 20. T 4 : Health state estimation and forecasting of ca- 
pacitance loss (%) at different times t p during the aging time; 
t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 


Aging time (hr) 

Figure 21. T 4 : Detail of the health state estimation and fore- 
casting of capacitance loss (%) at different times t p during the 
aging time; t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 
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* Measured o Filtered Predicted 


^ Measured o Filtered Predicted 



Figure 22. X 5 : Health state estimation and forecasting of ca- 
pacitance loss (%) at different times t p during the aging time; 
t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 


Aging time (hr) 

Figure 23. X 5 : Detail of the health state estimation and fore- 
casting of capacitance loss (%) at different times t p during the 
aging time; t p = [0, 24, 47, 71, 94, 116, 139, 149, 161, 171]. 
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A Prognostics Alpha-Lambda Performance 
Metric 




Figure 26. X 5 : Alpha-Lambda Prognostics Metric (A = 0.5 
and a = 0.3). 


Figure 24. T 2 : Alpha-Lambda Prognostics Metric (A = 0.5 
and a = 0.3). 



Figure 25. T 3 : Alpha-Lambda Prognostics Metric (A = 0.5 
and a = 0.3). 



Figure 27. T 5 : Alpha-Lambda Prognostics Metric (A = 0.5 
and a = 0.3). 
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Motivation 


• Electronic components have increasingly critical role in 
on-board, autonomous functions for 

- Vehicle controls, communications, navigation, radar systems 

• Future aircraft systems will rely more on electronic 
components 

• Assumption of new functionality increases number of 
electronics faults with perhaps unanticipated fault modes 

i a # i i i 1 r Line Replaceable Unit 

• We need understanding of m.pwrawi, 


behavior of deteriorated 
components to develop 
capability to anticipate 
failures/predict remaining RUL 





Images courtesy : Boeing 



2012 EPHM 



Prognostics Research Approach 

for Electronics 


iV 


2012 EPHM 



Model-based prognostics (1/3) 


x(t) = f(x(t),u(t )) + w(t) 
y(t) = h(x(t)),u(t)) + v(k) 


R(tp) — tEOL — t 


p 





RUL{t p ) 


State vector includes 
dynamics of the 
degradation process 

It might include nominal 
operation dynamics 

EOL defined at time in 
which performance 
variable cross failure 
threshold 

Failure threshold could 
be crisp or also a 
random variable 
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Model-based prognostics (2/3) 


Tracking of health „ 
state based on 
measurements 

Forecasting of 
health state until 
failure threshold is 
crossed 

Compute RUL as 
function of EOL 
defined at time 
failure threshold is 
crossed 


Measured 
O Filtered 
Predicted 
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Aging time (hr) 
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Fault Dimension (a) 


Model-based prognostics (3/3) 
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Methodology 


Test 

Trajectory 

%k = Axk-i + Bui -— i + w k — i 


y k = Hx k + 


{y(t 0 ), . . . ,y(t p )} 



RUL(t p ) 

— -.V 



Accelerated Aging Systems for 

Prognostics 
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Electrical Overstress Aging System 

Square wave amplified 
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Methodology 




Reduced 
order model 
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Dynamic modeling of the 
degradation process 
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Empirical degradation model 

• Based on observed degradation from 
capacitance parameter 

• Using training capacitor data to estimate 
degradation model parameters 

• Assumed exponential model based on 
capacitance loss 

• Parameter estimation with least-squared 
regression 

nt 

C k = e ah + P 

0 SO 100 190 200 

Aging Tim* <hr) 
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Degradation model results 



Validation 

test 

Test 

capacitor 

Training 

capacitor 

a 

(95% Cl) 

/3 

(95% Cl) 

9 

t 2 

#2 

#1, #3-#6 

0.0162 

(0.0160, 0.0164) 

- 0.8398 

(-1.1373,-0.5423) 

1.8778 

t 3 

#3 

#1, #2, #4 #6 

0.0162 

(0.0160, 0.0164) 

- 0.8287 

(-1.1211,-0.5363) 

1.9654 

Ti 

#4 

#l-#3, #5, #6 

0.0161 

(0.0159, 0.0162) 

- 0.8217 

(-1.1125, -0.5308) 

1.8860 

n 

#5 

#l-#4, #6 

0.0162 

(0.0161,0.0164) 

- 0.7847 

(-1.1134, -0.4560) 

2.1041 

n 

#6 

#l-#5 

0.0169 

(0.0167, 0.0170) 

- 1.0049 

(-1.2646, -0.7453) 

2.9812 


• The optimal parameter presented along the 95% confidence interval. 

• The residuals are modeled as a normally distributed random variable with zero 
mean and variance 
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Prognostics method and results 
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Prognostics algorithm 

• Implementation of prognostics algorithm with 
Kalman filter 

• Capacitance loss considered as state variable 

• EIS measurements and lumped parameter 
model used to obtained measured capacitance 
loss values 

• Empirical degradation model used to generate 
the state transition equation 

• Use one Capacitor for testing and the rest for 
model parameter estimation (leave on out test) 

• Failure threshold of 20% drop on capacitance 
based on MIL-C-62F 
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Kalman filter implementation 


State transition 
equation derived from 
degradation model 

at 

C k =e ah + (5 


\7 


dC 


dt 

Ct — Ct - At 

At 


= aC - a/3 


= aCt -At- aft 


Ct = (1 + (xA )C i - At — af) A 
Ck = (1 + aA. )Ck - 1 — cxfiAk 


State-space model for 
filter implementation 


Ck — AkCk - 1 + BkU + V 
yk = hCk + w, where 
Ak = (1 + A/), 

Bk = -a/3 A k , 
h = 1, u = 1. 
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Prediction mode 

• Assumed measurements are not available at 
some point in time 

• Filter used in forecasting mode to predict 
future states 

• Predictions done at 1 hr. intervals 

• State transition equation used to propagate 
state (n: number of prediction steps, /: last 
measurement at t t ) 

C, + „=A"C, + ^A'B 

i-0 
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Tracking and forecasting (Cap. #6) 
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Relative Accuracy 




RA = 100 1 - 


RUL* - RUL 


RUL 


* 


tp 

RAt2 

RA T3 

RAji 4 

RAt 5 

RAtq 

RA 

24 

94.8 

95.5 

91.9 

96.9 

99.7 

95.5 

47 

97.4 

99.3 

96.4 

96.7 

91.7 

96.7 

71 

87.5 

91.9 

84.5 

94.1 

97.1 

91.9 

94 

85.6 

90 

78.9 

94.8 

94.2 

90 

116 

86 

99.1 

76.5 

98 

96.2 

96.2 

139 

77.8 

95.8 

53.1 

96.7 

81.1 

81.1 

149 

82.1 

98.4 

46.9 

94.8 

86.6 

86.6 

161 

77.2 

87.3 

16.6 

87.5 

89.8 

87.3 

| 171 

26.6 

26.4 

N/A 

34.8 

63.7 

30.7 | 
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Contribution 


• RUL prediction algorithm based on accelerated 
life test data and an empirical degradation 
model 

- Identification of the lumped-parameter model for a 
real capacitor as a viab e reduced-order model for 
prognostics-algorithm development 

- Identification of ESR and C as precursor of failure 
feature parameters 

- Development of an empirical degradation model 
based on accelerated life test data which accounts 
for shifts in capacitance as a function of time 

- Implementation of a Bayesian based health state 
tracking and RUL prediction algorithm based on the 
Kalman filtering framework 
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Comments 

• Proposed approach requires further 
development 

- Results presented over accelerated aging test time 
scale 

- The empirical degradation model can be improved. 

- Degradation model and prediction algorithm assume 
constant loading and environmental conditions 

- Explore more sophisticated Bayesian tracking 
algorithms if required to handle variable loading and 
operational conditions as well as degradation models 
with time varying parameters 

- Uncertainty representation in the forecasting section 
and model uncertainty assessment under the 
Bayesian tracking framework 
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Thank You 

This work was funded by NASA, Aviation Safety 
Program, IVHM and SSAT Projects 
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Results (all) 
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RUL Results for all test cases 



tp 

RUL* 

RU Lrjp 2 

RU Lrp 3 

RU Lrp 4 

RU Lrjp^ 

RU Lrp 0 

24 

151.04 

158.84 

164.88 

158.76 

167.76 

159.89 

47 

128.04 

131.32 

134.08 

128.35 

135.32 

125.91 

71 

104.04 

117.01 

119.88 

115.37 

122.63 

116.41 

94 

81.04 

92.69 

96.64 

93.09 

97.6 

95.42 

116 

59.04 

67.28 

65.39 

67.77 

69.5 

65.71 

139 

36.04 

44.01 

44.72 

46.88 

49.4 

53.75 

149 

26.04 

30.67 

32.41 

33.55 

35.92 

39.95 

161 

14.04 

17.23 

18.28 

18.2 

22.64 

25.6 

171 

4.04 

1.07 

2.89 

N/A 

5.52 

8.45 
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Degradation on lumped parameter 

model 



C and ESR are estimated from 
EIS measurements 



Aging Time (hr) 
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